Radio-Frequency Spectroscopy of Ultracold Fermions 
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Radio-frequency techniques were used to study ultracold fermions. We observed the absence of 
mean-field "clock" shifts, the dominant source of systematic error in current atomic clocks based on 
bosonic atoms. This is a direct consequence of fermionic antisymmetry. Resonance shifts propor- 
tional to interaction strengths were observed in a three-level system. However, in the strongly inter- 
acting regime, these shifts became very small, reflecting the quantum unitarity limit and many-body 
effects. This insight into an interacting Fermi gas is relevant for the quest to observe superfluidity 
in this system. 
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Radio-frequency (RF) spectroscopy of ultracold atoms 
provides the standard of time. However, the reso- 
nance frequencies are sensitive to the interactions be- 
tween atoms, leading to the so-called clock shifts of the 
unperturbed resonances [1]. These shifts limit the accu- 
racy of current atomic clocks [2,3], but can also be used 
to characterize atomic interactions. 

RF spectroscopy has previously been applied to cold 
atoms to determine the size and temperature of atom 
clouds [4,5]. RF methods have also been used for evap- 
orative cooling, for preparing spinor Bose-Einstein con- 
densates (BEG) [6,7], and as an output coupler for atom 
lasers [5,8]. In all these experiments, shifts and broad- 
enings due to atomic interactions were negligible. Re- 
cently, density-dependent frequency shifts of RF tran- 
sitions were observed in rubidium [9] and sodium [10] 
BEGs. These frequency shifts are proportional to the 
difference in mean-field energies of two internal atomic 
states and allow scattering lengths to be extracted. Mean 
field shifts in BEGs have been observed also by optical 
spectroscopy [11,12]. 

Here, we apply RF spectroscopy to ultracold clouds of 
fermions and demonstrate several phenomena: (1) the 
absence of a clock shift in a two-level system due to 
fermionic antisymmetry, (2) the emergence of mean-field 
shifts in a three-level system after the relaxation of pair 
correlations, (3) the limitation of mean-field shifts due 
to the unitarity limit, and (4) the universality of the in- 
teraction energy in a dense cloud, independent of the 
attractive or repulsive nature of the two-particle interac- 
tions. 

Research in ultracold fermions has advanced rapidly, 
with six groups now having reached quantum degeneracy 
[13-18]. A major goal is to induce strong interactions by 
tuning magnetic fields to scattering resonances (called 
Feshbach resonances). Under these conditions, Gooper 
pairs of fermions may form, leading to superfluidity. This 
would establish a model system for studying Bardeen- 
Gooper-Schrieffer (BGS) pairing at densities nine orders 
of magnitude lower than in previous realizations in ^He 



and superconductors. We show that RF spectroscopy can 
be used to characterize interactions between fermions in 
the regime where superfluidity has been predicted [19,20]. 

Our experimental technique for preparing ultracold 
fermions has been considerably improved since our ear- 
lier work [17,21]. As the Pauli exclusion principle sup- 
presses elastic collisions between identical fermions at 
low temperatures and prevents evaporative cooling, we 
cooled fermionic ^Li sympathetically with bosonic ^^Na 
loaded into the same magnetic trap. In contrast to pre- 
vious work, we cooled both species in their upper hy- 
perfine states (^SNa : \F,mF) = |2,-l-2), ^Li : \F,mF) = 
]3/2, -f 3/2)). This led to a reduction of inelastic loss 
processes, boosting our final fermion atom numbers by 2 
orders of magnitude. We could produce BEGs contain- 
ing up to 10 million sodium atoms in the ]2, +2) state by 
evaporatively cooling pure bosonic samples in the mag- 
netic trap. For a Bose-Fermi mixture, the finite heat 
capacity of the bosons limited the final lithium temper- 
ature after the 30 s evaporation cycle to ^0.3rF for 10 
million fermions and Tp for 50 million fermions [22], 
where Tp is the Fermi temperature. 

The spin states of ®Li of most interest for superfluid 
pairing are the two lowest states ]1) and ]2) (]l/2, -t-1/2) 
and [1/2, —1/2) at low field), which are predicted to have 
an inter-state s-wave Feshbach resonance at ~ 800 G 
[23,24]. However, both states are high-field seeking at 
these fields, making them unsuitable for magnetic trap- 
ping. We therefore transferred the atoms into an opti- 
cal trap. For these experiments, 6-8 million ]3/2, -1-3/2) 
lithium atoms were loaded into the optical trap at T 
Tp ~ 35 /xK [25] . The atoms were then transferred to 
the lowest energy state ]1), using an adiabatic frequency 
sweep around the lithium hyperfine splitting of 228 MHz. 
DG magnetic fields of up to ~ 900 G were applied, a 
range encompassing the ] 1) — ]2) Feshbach resonance. Us- 
ing RF-induced transitions near 80 MHz, we could create 
mixtures of states Jl), [2) and [3) (13/2,-3/2) at low 
field) , and explore interactions between fermions in these 
states. 
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Collisions between atoms cause a shift of their energy, 
which is usually described by the mean-field effect of all 
the other atoms on the atom of interest. For example, 
atoms in state |2) experience an energy shift "^^^^ »-iai2 
due to the presence of atoms in state Here h is 

Planck's constant divided by 27r, m is the mass of the 
atom, m is the density of |1) atoms and ai2 is the inter- 
state scattering length. We use the convention that pos- 
itive scattering length corresponds to a repulsive inter- 
action. Density-dependent shifts of the resonance fre- 
quency for the transition connecting two states have been 
observed in laser-cooled [1] and Bose-condensed clouds 
[9,10]. 

In the case of ultracold fermions, only interactions be- 
tween different internal states are allowed. For a system 
of density n, let us compare the energy of a gas pre- 
pared purely in state to a gas in which one atom 
is transferred into state |2). The energy difference is 
hi'12 + '^^jH na 12, where is the rc^sonancc; frequency 
of the non-interacting system. Similarly, the energy dif- 
ference between a gas prepared purely in state |2), and 
a gas in which one atom is transferred into state |1) is 

hi'12 - ^nai2. 

However, these energy shifts should not affect the res- 
onance for a coherent transfer out of a pure state. For 
fermions in the initial pure state, the pair-correlation 
function vanishes at zero distance due to the antisym- 
metry of the wavefimction. During any coherent trans- 
fer process, the state vectors of all the atoms rotate "in 
parallel" in Hilbert space, i.e. the superposition of the 
two spin states has the same relative phase for all atoms. 
Thus, the atoms remain identical and cannot interact in 
the s-wave regime. The mean-field energy is thus estab- 
lished only after the coherence of the; superposition state 
is lost and the pair correlations have relaxed, forming a 
purely statistical mixture of the two states. 

It is a consequence of Fermi statistics that spectro- 
scopic methods do not measure the equilibrium energy 
difference between the initial and final state of the sys- 
tem, but rather the unperturbed resonance frequency. 
The expected absence of the clock shift has led to sug- 
gestions for the use of fermions in future atomic clocks 
[26]. Our work presents an experimental demonstration 
of this phenomenon. 

We determined the transition frequency between states 
|1) and 1 2), first starting with a pure state and then 
with a pure state |2) sample. The absence of a splitting 
between these two lines proves the suppression of the 
clock shift. Fig. 1 shows an example of such measure- 
ments. The magnetic field was ramped up to 570 G with 
the cloud in state |1). At this field, a\2 ~ 150 ao- There- 
fore, the expected equilibrium mean-field shifts were 
Av = ±5 kHz for our mean density of 3 x lO^^ ,,^-3 [27] . 
The interaction between states |1) and |2) at this mag- 
netic field was also observed in the mutual evaporative 



cooling of the two states in the optical trap. RF pulses 
of 140 fis duration were applied at frequencies near the 
unperturbed resonance z/i2 ~ 76 MHz. Atoms in states 
|1) and 1 2) could be monitored separately by absorption 
imaging, since they are optically resolved at this field. 
We observed a suppression of the clock shift by a factor 
of 30 (Fig. 1). Using the same method, the absence of 
the clock shift was observed at several other magnetic 
fields. In particular, we observed a suppression of more 
than three orders of magnitude at 860 G [28] . 
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FIG. 1. Absence of the clock shift. RF transitions were 
driven between states |1) and |2) on a system prepared purely 
in state |1) (filled circles), and purely in state |2) (open cir- 
cles). Mean-field interactions would result in 5 kHz shifts for 
the two curves in opposite directions. Gaussian fits (solid 
lines) to the data are separated by 0.04 ±0.35 kHz. This gives 
a clock shift suppression factor of 30. 

P-wavc interactions could lead to a non- vanishing clock 
shift. However, at these low temperatures they are pro- 
portional to T otTf, whichever is higher, and are there- 
fore strongly suppressed. 

Mean-field shifts and scattering lengths can be ob- 
served spectroscopically by driving transitions from a 
statistical mixture of two states to a third energy level. 
(While this work was in progress, use of a similar method 
to measure scattering lengths in fermionic was re- 
ported [29].) Specifically, wc recorded the difference be- 
tween the RF spectra for the |2) |3) transition in the 
presence and in the absence of state |1) atoms. The pres- 
ence of atoms in state |1) is then expected to shift the 
resonance by [30]: 

Ai'= — ni(ai3-ai2) (1) 
m 

In our experimental scheme to determine the interac- 
tion energy at different magnetic fields (Fig. 2), the sys- 
tem was prepared by ramping up the magnetic field to 
500 G with the atoms in state |1). Either partial or com- 
plete RF transfer to state |2) was then performed. The 
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number of atoms in state |1) was controlled by adjusting 
the speed of a frequency sweep around the |1) ^ |2) res- 
onance. A fast, non-adiabatic sweep created a superpo- 
sition of the two states, whereas a slow, adiabatic sweep 
prepared the sample purely in state |2). A wait time of 
200 ms (see below) was allowed for the coherence between 
states |1) and |2) to decay and the system to equilibrate. 
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FIG. 2. Schematic of the mean-field measurement, (a) Hy- 
perfine structure of the ground state of ''Li. (b and c) Ex- 
perimental scheme: (b) preparation of a mixture of atoms in 
states |1) and |2), and (c) RF spectroscopy of the |2) — » |3) 
transition. 

Typical parameters for the decohered |1) — |2) mixture 
were mean-density n\ ~ 2.4 x 10^^ cm~^ and T ^0.7 Tp- 
The magnetic field was then changed to the desired value, 
and the transition from state |2) to state |3) was driven 
with 140 /xs RF pulses (Fig. 2(c)). We monitored the ap- 
pearance of atoms in state |3) and the disappearance 
of atoms from state |2), using simultaneous absorption 
imaging. Fig. 3(a) shows the unperturbed and perturbed 
resonances at B = 480 G. The position of the unper- 
turbed resonance i/23 also determines the magnetic field 
to < 0.1 G accuracy. Fig. 3(b) shows absorption images 
of atoms in state |3), obtained for different values of the 
applied radio-frequency. One can clearly see the spatial, 
and thus the density dependence of the mean-field shift: 
close to the unperturbed resonance, the low density wings 
of the cloud are predominantly transferred, whereas the 
high density central part of the cloud is transferred only 
at sufficient detuning. To suppress spurious effects from 
this spatial dependence, only a small central part of the 
images was used to extract the transferred atomic frac- 
tion. 
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FIG. 3. Representative spectra at 480 G. (a) Fraction of 
atoms transferred from |2) to |3), with |1) atoms absent (filled 
circles), and present (open circles). The mean-field shift is 
computed from gaussian fits to the data (solid lines), (b) 
Spatial images of state |3) for the perturbed resonance. The 
optical trap was turned off immediately after the RF pulse 
and absorption images of the atoms were taken after 120 /its 
expansion time. The central section of ~ 150 /itm vertical 
extent was used to extract the transferred fractions in (a), (b) 
also shows images of states |2) and |1) for zero RF detuning. 
States |3) and |2) were imaged simultaneously to observe their 
complementary spatial structure. State |1) was imaged after 
760 /is expansion time to record its density for normalization 
purposes. 

To ensure that our mean-field measurements were 
performed on a statistical mixture, we measured the 
timescale for decoherence in our system. The decay of 
the |1) — 1 2) coherence at 500 G was observed by mon- 
itoring the |2) — > |3) transfer at the measured unper- 
turbed resonance 1^23, as a function of wait time (Fig. 4). 
For wait times small compared to the decoherence time 
of the |1) - |2) superposition, the |2) ^ |3) RF drive 
places each atom in an identical three-state superposi- 
tion. All mean-field shifts are then absent and the result- 
ing transfer is unchanged from the unperturbed case. For 
longer wait times, the |1) — |2) superposition decoheres 
and mean-field interactions set in. This shifts the reso- 
nance frequency of the |2) |3) transition, reducing the 
transferred fraction at 1^23- The measured decoherence 
time of ~ 12 ms was attributed mainly to the sensitivity 
of z/12 to magnetic field variations across the cloud. These 
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inhomogencitics cause the relative phase of the |1) — |2) 
superposition in different parts of the trap to evolve at 
different rates, given by the local ui2- Atoms which travel 
along different paths within the trap therefore acquire 
different phases between their |1) and |2) components. 
Being no longer in identical states, s-wave interactions 
between them are allowed. The inhomogencitics scale 
with the applied magnetic field B while the sensitivity 
of the transition scales with dv\2/dB. We would thus 
expect the dccohcrence time to vary inversely with the 
product of these two quantities. Our hypothesis is sup- 
ported by our observation of longer decoherence times at 
higher fields, where B x dvyijdB is lower. 
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FIG. 4. Emergence of mean-field shifts due to decoherence 
at 500 G. Decoherence leads to a reduction of the |2) — > |3) 
transfer at the unperturbed resonance 1*23. An exponential fit 
to the data (solid line) gives a time constant of 12 ms. 

Fig. 5 summarizes the results of our mean-field mea- 
surements for a wide range of magnetic fields up to 750 G. 
For magnetic fields up to 630 G, our data can be ex- 
plained fairly well using Eq. 1 with the theoretical calcu- 
lations of the scattering lengths shown in Fig. 6, and an 
effective density of n\ = 2.2 x lO^^cm"'^, consistent with 
the initial preparation of the system at 500 G. A narrow 
resonance of ai2 at ~ 550 G [21,24,31] is indicated by the 
data, but was not fully resolved. We also see additional 
structure near 470 G, which is not predicted by theory 
and deserves further study. 



20- 



TTlO- 



-20- 



I 






/ % 























400 500 600 

Magnetic Field [G] 



700 



FIG. 5. Frequency shift vs. magnetic field for the |2) — > |3) 
resonance due to atoms in state |1}. The shifts are computed 
by monitoring the arrival fraction in state |3) for 140 /is RF 
pulses, except at 750 G. Here, because of strong inelastic losses 
between |3) and |1) atoms, we monitored the loss of atoms in 
state 1 2) after applying RF sweeps of 3 ms duration and 2 kHz 
width. All the data points are normalized to the same atom 
number in state jl). The fit at low fields (solid line) uses Eq. 1 
with m = 2.2 X 10^"^ cm""^ and the theoretical calculations of 
the scattering lengths. The error bars reflect uncertainty in 
the state |1) atom number, and the uncertainty in the gaus- 
sian fits to the spectra. The dashed line indicates the position 
of the predicted 013 resonance. 

The most important results of this paper are our obser- 
vations for fields above 630 G, where the measured shifts 
strongly deviate from the predictions of Eq. 1. In the re- 
gion between 630 G and 680 G, the two scattering lengths 
are expected to be large and positive, with ai3 ^ a\i 
(Fig. 6). Eq. 1 would thus predict large positive mean- 
field shifts. In contrast, we observe very small shifts, 
indicating almost perfect cancellation of the two contri- 
butions. Even more surprisingly, we observe essentially 
no mean-field shifts between 680 G and 750 G, where the 
two scattering lengths are predicted to be very large in 
magnitude and of opposite signs, and in a simple picture 
should add up to a huge negative shift. These results are 
evidence for new phenomena in a strongly interacting 
system, where the scattering length becomes comparable 
to either the inverse wavevector of interacting particles, 
or the interatomic separation. 

Eq. 1 is valid only for low energies and weak interac- 
tions, where the relative momentum of the two particles, 
satisfies -C f/|fl|. For arbitrary values of fca, the 
s-wave interaction between two atoms is described by re- 
placing the scattering length a with the complex scatter- 
ing amplitude: 



/ = 



1 -h A:2a2 



(1 — ika) 



(2) 



4 



The real part of / determines energy shifts, and hence 
the ground state properties of an interacting many-body 
system. The imaginary part determines the (inverse) 
Ufctime for clastic scattering out of a momentum state, 
and hence the dynamic properties of the system such as 
thermalization rates. For k\a\ — > oo, the elastic cross- 
section (T = 47rIm(/)/fc monotonically approaches the 
well known "unitarity limited" value of Air/k'^. On the 
other hand, the two particle contribution to the mean- 
field energy, proportional to — Rc(/) = a/(l + fc^a^), 
peaks at \a\ = 1/k, and then, counter-intuitively, de- 
creases as l/\a\ for increasing \a\. It has been shown [32] 
that averaging Re(/) over a zero-temperature Fermi dis- 
tribution with Fermi momentum nkp, limits its absolute 
value to l.Ob/kp, and drastically weakens its dependence 
on the exact value of a in the /cpjaj > 1 regime. This re- 
sults in a prediction for the mean-field energy which is 
sensitive to the sign of the scattering length, never ex- 
ceeds OA^Ep, and remains finite for kplal ^ 1- Hence, 
this approach could qualitatively explain our results in 
the 630 — 680 G region, but it is in clear contradiction 
with negligible resonance shifts in the 680 — 750 G region 
[33]. 
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FIG. 6. s-wave scattering lengths ai2 and ais as a function 
of magnetic field, obtained from a highly model-independent 
quantum scattering calculation. The calculation makes use of 
the presently available ®Li experimental data [39] in a coupled 
channel approach to deduce accumulated phases that charac- 
terize the less well-known short range parts of the ^Li -|- ®Li 
scattering potential [31]. ai2 has a narrow Feshbach resonance 
at 550 G and a wide one at 810 G. ais has a wide Feshbach 
resonance at 680 G. 

We suggest that these discrepancies might be due to 
the fact that we arc in the high density regime, where 
n|ap approaches unity. In a degenerate Fermi gas, the 
interparticle spacing is comparable to the inverse Fermi 



wavcvcctor, kp = Ott^tt,. Hence, the unitarity limit co- 
incides with the breakdown of the low-density approxi- 
mation {n\a\^ <C 1) and higher-order many-body effects 
can become important. Some recent many-body calcu- 
lations [34-36] suggest that in the regime kpla] ^ 1 (or 
n|a|^ » 1), the interaction energy is always negative and 
independent of both sign and magnitude of a. This sug- 
gests that whenever the scattering length is large, either 
positive or negative, the interaction energy is a univer- 
sal fraction of the Fermi energy [32]. This is a possible 
explanation for the small line shifts that we observed for 
fields above 630 G, where the interactions are strong in 
both states. 

This picture is consistent with other recent experimen- 
tal observations [29,32,37,38]. Expansion energy mea- 
surements in a mixture of states |1) and |2) of ^Li [38], 
showed a negative interaction energy at 720 G, which is 
on the repulsive side of the predicted Feshbach resonance. 
RF spectroscopy in ^"K [29] has also shown some sat- 
uration of the mean-field in the vicinity of a Feshbach 
resonance, which may reflect the unitarity limit. 

In characterizing an interacting Fermi gas by RF spec- 
troscopy, we have demonstrated absence of clock shifts in 
a two-level system, and introduced a three-level method 
for measuring mean-field shifts. For strong interactions, 
we have found only small line shifts which may refiect 
both the unitarity limit of binary collisions and many- 
body effects. It would be very important to distinguish 
between two-body and many-body effects by studying 
the gas over a broad range of temperatures and den- 
sities. In a very dilute and very cold gas, the weakly 
interacting limit could be extended to very large values 
of |a|, thus allowing for direct verification of molecular 
calculations. This presents experimental challenges, be- 
cause cooling changes the density and the temperature 
together. It would also be interesting to study similar 
phenomena in bosonic gases, in order to distinguish to 
what extent the high density many-body effects depend 
on quantum statistics. This new insight into the physics 
of strongly interacting Fermi gases must be taken into 
account in the search for superfiuidity in these systems. 
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